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Abstract: We review some recent results concerning the quantitative analysis of the universality classes 
of two-dimensional statistical models near their critical point. We also discuss the exact calculation of the 
two-point correlation functions of disorder operators in a free theory of complex bosonic and fermionic 
field, correlators ruled by a Painleve differential equation. 



1. Introduction 

One of the most important successes of Quantum 
Field Theory (QFT) in recent years is the quanti- 
tative analysis of the universality classes of two- 
dimensional statistical mechanical models near their 
second order phase transition points. This has 
been possible, in particular, thanks to new devel- 
opments in the methods of computing off-critical 
correlation functions. The aim of this talk is to 
review some of the recent results relative to these 
two important subjects. Our discussion is divided 
in two parts: in the first part, we shall see how 
correlation functions can be used to extract useful 
numbers (the so-called universal ratios) that can 
be directly compared with experiments. In the 
second part, we shall concentrate on a particularly 
interesting class of models, those of fermionic and 
bosonic theories with a symmetry, for which 
there is an elegant approach to derive differential 
equations satisfied by the two-point functions of 
their disorder operators. 

2. Universal Amplitude Ratios 

Consider a statistical model with n relevant fields 
Lpi{x) at criticality. Near the critical point, its 
action can be parameterised as 



CFT 



+ 9i I ^i{x)d^x 



(2.1) 



where Ac ft is the action which encodes the data 
of the Conformal Field Theory of the critical point. 

*Talk presented by G. Mussardo. 



The short -distance behaviour of the two-point func- 
tions is given by 



1 



|4Ai 



for I a; I — s- 0, where is the conformal dimen- 
sion of the operator tpi. Hence, their conjugate 
couphng constants gi behave as gi ~ p^D~2Ai ^ 
where A is a mass scale. Away from criticality 
there will be generally a finite correlation length 
f which in the thermodynamical limit scales as 
^ ~ a {Kigi) , where a ~ A~^ may be 

regarded as a microscopic length scale. The di- 
mensionless quantities are non-universal met- 
ric factors specific of the system representing the 
universality class and depending on the units cho- 
sen for measuring the external sources gi. If there 
are several deformations of the conformal action, 
the most general expression for the scaling form 
of the correlation length may be written as 



e = 6 = a(if,.9,)-™T£ 



K 



j9j 



where 



_D-2A, 



{K,gi)^^ 

(2.2) 

are the so-called crossover 



exponents whereas Ci are universal homogeneous 



scaling functions of the ratios 



There are 



of course several (but equivalent) ways of express- 
ing these scaling forms, depending on which cou- 
pling constant is selected as a prefactor. In the 
limit where gi ^ Q {I ^ i) but 7^ 0, equation 
(|2.2(l becomes 



Si 



(2.3) 
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Let us consider now the free-energy . . . , £f„], 
a dimensionless quantity defined by 

e-/(si.-..,sn) =y D^e-[-4c^'r+Er=iS'/'^*(^^'^°^ 

(2.4) 

Under the hypothesis of hypcrscahng, its singular 
part (per unit of volume) f[gi, . . . ,gn] is propor- 
tional, in the thermodynamical limit, to the D-th 
power of the correlation length. Depending on the 
scaling form adopted for the correlation length, 
we have correspondingly several (but equivalent) 
ways of parameterizing this quantity 



fi[gi,...,gn] = {Kigi)T^ Ti 



{Kigi)^ 



(2.5) 

The functions are universal homogeneous scal- 



ar 



j9j 



ins functions of the ratios , , . . 

Prom the above expression we can derive the 
expectation values of several quantities. We will 
use the notation (...), to denote expectation values 
computed in the off-critical theory obtained by 
keeping (at the end) only the coupling constant 
gi different from zero. Basic quantities are the 
vacuum expectation values (VEV) of the fields ipj 
which can be parameterized as 



/ \ 



91=0 



with 



C-2Ai 



(2.6) 



(2.7) 



The above relations can be equivalently expressed 
as 

9i = Dii{{iPj)i)^ , 



with 



1 



(2.8) 
(2.9) 



The generalized susceptibilities of the model are 
defined by 



jk 



d 

dgk 



dgkdgj 



(2.10) 



S!=0 



They are obviously symmetrical in the two lower 
indices. By extracting their dependence on the 
coupling constant gi, they can be expressed as 



'-jk 



^k9i 



(2.11) 



with 



2Aj-+2Afc-0 



(2.12) 



The various quantities obtained by taking the deriva- 
tives of the free-energy obviously contain the quan- 
tities Ki which make their values not universal. 
However, it is easy and always possible to con- 
sider special combinations thereof in such a way 
to cancel out all metric factors. Universal Ratios 
usually considered in the literature are those given 
below 



ii jk 




(2.13) 

(2.14) 
(2.15) 

(2.16) 
(2.17) 



From their definition, these quantities are pure 
numbers attached to the universality classes and 
therefore they can be used to characterize them. 
Contrary to the critical exponents (which are char- 
acteristic of the critical point) they carry informa- 
tion about the scaling region. Moreover, the am- 
plitude ratios are numbers which typically present 
significant variations between different classes of 
universality, whereas the critical exponents usu- 
ally assume small values which only vary by a 
small percent when changing the universality class. 
Hence the universal ratios may be ideal marks 
of the critical scaling regime []. It is also worth 
emphasizing that, from an experimental point of 
view, it should be simpler to measure universal 
amplitude ratios rather than critical exponents: to 
determine the former quantities one needs to per- 
form several measiircments at a single, fixed value 
of the coupling which drives the system away from 
criticality whereas to determine the latter, one 
needs to make measurements over several decades 
along the axes of the off critical couplings. More- 
over, although not all of them are independent, 
the universal ratios are a larger set of numbers 
than the critical exponents and therefore permit 
a more precise determination of the class of uni- 
versality. 
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3. Quantum Field Theory Approach 

It is interesting to see how in two dimensions Quan- 
tum Field Theory provides a powerful approach 
to compute the Universal Ratios of a given class 
of universality in a non-perturbative way. In the 
following we concentrate mostly our attention on 
the off-critical models obtained as deformations of 
the so-called Minimal Models of two-dimensional 
Conformal Field Theory [] . 

As we have seen before, each coupling con- 
stant Qi relative to the relevant operator tfi{x) 
of the model under investigation is a dimensional 
quantity which can be related to the lowest mass- 
gap rrii — of the off-critical theory according 
to the formula 



rrii = d g- 



(3.1) 



When the QFT associated to the action (|2.1|l is in- 
tegrable, the pure number Ci can be exactly deter- 
mined by means of the Thermodynamical Bethe 
Ansatz []. When the theory is not integrable, the 
constant Ci can be nevertheless determined by a 
numerical method, based on the so-called Trun- 
cated Conformal Space Approach []. In conclu- 
sion, for all individual deformations of a given 
model, we are able to completely set the rela- 
tionship which links the coupling constant to the 
mass-gap of the theory and therefore to switch 
freely between these two variables. 

Additional quantities which can be determined 
by QFT are the matrix elements of the order pa- 
rameters, the simplest ones being the vacuum ex- 
pectation values (VEV). They are parameterised 



as 



(3.2) 



When the theory is integrable, the constant Bji 
can be fixed exactly, thanks to the results of a 
remarkable series of papers []. When it is not 
integrable, the constant Bji can be nevertheless 
estimated by means of a numerical approach, as 
firstly shown in []. Hence, we are able to deter- 
mine completely also these quantities. Moreover, 
as shown in [] , a generalization of the numerical 
approach of ref. [] often leads to a reasonable esti- 
mate of the matrix elements of the order param- 
eters between the vacuum states and some of the 
excited states, as for instance {Q\ipj\Ai^)i where Ak 



is a one-particle state of mass Mk. These quan- 
tities turn out to be useful for obtaining sensible 
approximation of the large-distance behavior of 
several correlators. 

Another important and useful piece of infor- 
mation on the off-critical models can be obtained 
by exploiting the properties of the stress-energy 
tensor T^i,{x). In the presence of the perturbing 
field ipi, the trace of the stress-energy tensor is 
given by 

Q{x) = 2715,(2 -2 A,) . (3.3) 
It enters a useful sum rule, the so called A-theorem 



1 



47r((y5 



d2x(e(x)^,(0))^ , (3.4) 



which relates the conformal dimension A^ of the 
field (fj to its VEV and to the integral of its con- 
nected off-critical correlator with Q{x). It is easy 
to see that the above formula simply expresses 
the content of the fluctuation-dissipation theo- 
rem and when the above integral diverges, so does 
the VEV in the denominator, in such a way that 
eq. (|3.4|) always keeps its validity []. 

As discussed above, basic quantities entering 
the universal ratios are the generalized suscepti- 
bilities r*j.. By using the fluctuation-dissipation 
theorem, they are given by 



(3.5) 



By extracting their dependence on the coupling 



constant gi, we have f^^^, — Tjj^g- 

1 



T]^ = 0^^^+^^"-^ I dr 



with 



(3.6) 

Some of the above susceptibilities can be deter- 
mined exactly, such as the components r*^, , whose 
values are provided by the A-theorem sum rule 

r^fc = — i^^ki ■ (3.7) 

J- - Afe 

In all other cases, when an exact formula is not 
available, the strategy to evaluate the generalized 
susceptibilities relies on two different possible rep- 
resentations of the correlators. As we will see, 
these representations have the advantage of con- 
verging very rapidly in two distinct regions and to 
smoothly join in between. 
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The first representation is based on Conformal 
Perturbation Theory (with the use of the non- 
analytic expression of the VEV) []. In this ap- 
proach the two-point correlators are expressed as 

{v^{x)vM) = 5]Cf,(ff;a;)(A,(0)) (3.8) 



where the structure functions Cji,{g;x) admit the 
expansion 



ri=0 



r 

(3.9) 

(r =1 a: |) and can be computed perturbatively 
in g, the coupling conjugated to the perturbing 
operator ip. For instance, the first order correc- 
tion to the conformal structure constant C' 



determined by [] 



p{0) 

hi 



C: 



'p(i) 



d^W {AP{(X,)ip{w)ip^{l)ip,{0))cFT , 

(3.10) 

where the prime indicates a suitable infrared (large 
distance) regularization of the integral. This rep- 
resentation allows a very efficient estimate of the 
correlation function in its short distance regime 
r < 

The second representation is based on the Form 
Factors (which will be discussed in more detail be- 
low) and allows an efficient control of the large dis- 
tance behavior, i.e. when r 3> ^. In this second 
representation, one makes use of the knowledge 
of the off-critical mass spectrum of the theory to 
express the correlators as 

oo 

(^.(x)^,(0))=^5l?M , (3.11) 



71=0 



where 



Jei>e2...>e 

X(0|(^,(0)| . ..AaM){- ■ ■ AaAen)\^jiO)\0) . 

I^ai (^i) • ■ • ^a„ {(^n)) ^rc the multi-particle states 
relative to the excitations of mass m^j, with rela- 
tivistic dispersion relations given by E — rrik cosh 0, 
p — mfeSinhS, where 9 is the rapidity variable. 
The spectral representation (|3.12l) is an expansion 
in the parameter e~ 5 , where = mi is the low- 
est mass-gap. 



Basic quantities of the large distance approach 
arc the Form Factors (FF), i.e. the matrix ele- 
ments of the operators fi on the physical asymp- 
totic states 



ai 



d9i dOn -rV" mfccosliffc 



27r 27r 



„) = {o\ipm\AaAei) ■ ■ ■ AaM) 

(3.12) 

For scalar operators, relativistic invariance requires 
that the FF only depend on the rapidity differ- 
ences 9i — 9j . In integrable quantum field theories 
the FF are solutions of some functional equations 
[] which will be recalled in the next section. It is 
worth to stress that for most of the operators, it 
is sufficient to insert into the spectral representa- 
tions only their one-particle and two-particle FF 
in order to reach a very reasonable estimate of the 
correlators. 

In conclusion, both the expansions (ESJ and 
H3.11(l converge very rapidly in their domain (see, 
for instance []) and therefore it is possible to es- 
timate the integral H3.5|l by following this proce- 
dure: 

1. Express the integral in polar coordinates as 

r+oo 

f},-2^ / drr{ip,{r)ipkmi , (3.13) 

and split the radial integral into two pieces 
as 

r+oo 

/ drr{ip,{r)MO)yc 
Jo 

rR p+oc 

drr{...)l + / drr{...)l 

10 JR 

= h{R) + h{R) . (3.14) 

2. Use the best available short -distance repre- 
sentation of the correlator to evaluate /i (i?) 
as well as the best available estimate of its 
large-distance representation to evaluate I2 {R) . 

3. Optimize the choice of the parameter R in 
such a way to obtain the best evaluation of 
the whole integral. In practice, this means 
looking at that value of R for which a plateau 
is obtained for the sum of Ii (R) and I2 (R) ■ 
Said in another way, R belongs to that in- 
terval where there is an overlap between the 
short-distance and the long-distance expan- 
sion of the correlator. 
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The above methods have been apphed to com- 
pute universal ratios of several interesting univer- 
sality classes, such as the one of the Ising model 
[] , the Tricritical Ising Model [] , the Self- A voiding 
Walks [], the q-state Potts model and percolation 
[], to quote few. These field theoretical results 
have been confirmed through direct lattice esti- 
mates (see Refs. []). 

4. Disorder Operator and Fermion- 
Boson Correspondence 

From the above discussion, we have seen that a 
great deal of information can be extracted for the 
scaling region of a model also in the absence of 
an exact expression of the correlation functions. 
It is however an important theoretical problem to 
see whether it is possible to obtain a closed and 
exact formula for these quantities. It goes with- 
out saying that, in general, this is quite a difficult 
problem, partially solved only for few models. In 
the second part of this talk, we shall discuss a par- 
ticularly interesting example of such a computa- 
tion, which leads to a differential equation ruling 
the correlators. Our presentation is based on the 
results presented in the paper []. 

In order to start our discussion, let us consider 
the (euclidean) action 

[ d^x[df,(j)*d^q^+m^<j)*^+Wiq^*,cb)], (4.1) 



where the potential W is invariant under the Z 



N 



transformation (j) — > e^"/^0, (/>* 



-2iTT 



Let |0j), j ^ 0,1, . . . , N — I, he the N vacuua of 
the broken phase of this theory, with fikix), k = 
1, . . . ,N — 1 , the disorder operators which create 
the excitations (kink states) interpolating between 
the vacua |0j) and \Oj-^-k{mod n)) ■ These opera- 
tors carry k units of topological charge and satisfy 
the conjugation relation /i^(a;) = ^N-k(x). The 
mutual non-locality between the order and disor- 
der operators gives rise to a phase factor when 
they go around each other in the euclidean plane 
(z = xi+ ix2, z = Xi — ix2) : 



(ze2*'^,ze-2*'^ 



0*(ze2 



e « (f){z,z) Hk{0,0) , 
e-"-^<t>*{z,z) Mfc(0,0), 



These relations are valid along any Z^v-invariant 
renormalization group trajectory flowing out of 



the phase transition point and characterise the op- 
erators /ifc {x) beyond their initial role of kink cre- 
ation operators in the broken phase. The theory 
H4.1|) with = is a particularly simple example 
of such a trajectory for which, as we will see, the 
correlation functions of disorder operators remain 
non-trivial due to the non-locality with respect to 
the boson. 

The correlation functions of non-local opera- 
tors in two-dimensional free massive theories were 
extensively studied from the point of view of the 
isomonodromy theory of differential equations in 
Ref. [] and a series of related papers. Several works 
(see for instance Refs.[]) have been devoted after- 
wards to dealing with this problem through more 
direct and general methods of quantum field the- 
ory. As we have seen above, an important ap- 
proach to compute the correlation functions is the 
one based on the form factors, in which the cor- 
relators are expressed as sums over multiparticle 
asymptotic states. In our example this series can 
be explicitly resummed and the correlators are 
given in terms of solutions of non-linear differen- 
tial equations. Moreover, there is a simple unified 
treatment of the bosonic and fermionic cases, with 
the final result expressed by the correspondence^ 

(A.(^)A^(0)> = ,,y , ,V ' (4-2) 

(V,/jv(a;)Vfc/Ar(0)) 

where the correlators on the l.h.s. are computed 
in the theory H4.1|l with W — 0, and those on the 
r.h.s. refer to the operators 14, (x) = exp[i-\/47r a (^(x 
in the sine-Gordon theory 

'1 



A 



sG 



d^X 



■ d^ifd'^if — /X cos l3ip 



(4.3) 



with P — \/At: and fi a suitably chosen mass scale. 
The theory (|4.3|) with /3 = y/An is in fact a free 
fermionic theory [] . For generic values of (3 the el- 
ementary excitations are the solitons and antisoli- 
tons interpolating between adjacent vacua of the 
periodic potential. Their interaction is attractive 
as long as (3 < y/An and they form topologically 
neutral bound states, the lightest one being the 
particle interpolated by the bosonic field in the 
action H4.3|) . These neutral particles are absent 
from the spectrum of asymptotic states in the re- 
pulsive region /3 > VAtt, and at the point (3 — \/4tt 
where the solitons behave as free Dirac fermions. 

^We use the notation 3>(a;) = $(a:)/(<J>). 
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Since the solitons are non-local in terms of the 
field ipix), the evaluation of both sides of Eq. H4.2() 
amounts to computing correlation functions of op- 
erators which are non-local with respect to non- 
interacting particles, bosons for the l.h.s. and 
fermions for the r.h.s. The different statistics of 
the two particles is responsible for the inversion in 
Eq. 

Consider then a theory of free, charge conju- 
gated particles A and A with mass m, and denote 
by ^a{x) a scalar operator having a non- locality 
phase e^''^" (e~^*^") with respect to (the field 
which interpolates) the particle A {A). Acting on 
the vacuum state |0), such an operator produces 
neutral states consisting of pairs A{6)A{(3). The 
equations satisfied by the form factors of these 
theories 

/,'^(0i,...,0„,/3i,...,/3„)= (4.4) 
= (O|$„(O)|A(0i), . . . , A(0„), A(/3i), . . . , A(/3„)) 

are a particular case of those holding for generic 
integrable theories (see e.g. []) 

/"(6'i, . . . ,6i, 6'i+i, . . . ,6n, Pi, ■ ■ ■ , (3n) = 
= S fn{di, . . . , 9i+i,6i, . . . ,6n, Pi, . . . , Pn), 

/;^(0i+2z7r,02,...,^„,/3i,...,/3„) = 
= S e^*'""/"(0i, . . . ,On,Pi,. . . ,(in), 
Resei_/3i=i^/"(6'i, . . . , . . . , /?„) — 

= tS^-Hl - e"'"")/„"_i(02, ..,e„,/32,. .,/?„), 

where 

1 for free bosons 
— 1 for free fermions . 



5" = 



(4.5) 



It is easy to check that the solution of the above 
system of equations is given by 

/„"(0i, . . . , . . . , /3„) = 5«("+2)/2(_ sinTra)" 

xe("-i^-'OEL,(^.-ft) , (4.6) 

where An is a n x n matrix (Aq = 1) with entries 

1 



A, 



cosh 



(4.7) 



and |^n|(5) denotes the permanent^ of An for S = 
1 and the determinant of An for S* = — 1. 

^The permanent of a matrix differs from the deter- 
minant by the omission of the alternating sign factors 
(-1) 



The two-point correlators are then given by 

1 



(n!)2 (27r) 



d01... dOndPl ...df3n 5^"'"'^(i \9l,...,Pn), 



where 

9i"'"'Ht)^f^idi,...,f3n)f:'iPn 

= (S-siuTra sin7ra')"e<"""')Sr=i(«'-ft) (4., 
l(S) 

and 



i = m|x| , e„ = ^(cosh^fc -I- cosh/3fe) 



Analogously to the procedure followed in [], we 
define a new n x n matrix M„ with entries 

M,j = M(e^,l3j) = (siuTra sinTra')^/^ 6"^'=°"''^' 



cosh 



(4.9) 



where 

and rewrite ) as 



Mn 

Mj 



(4.10) 



{S) 



i+i 



Finally we symmetrise with respect to the two sets 
of rapidities by introducing a charge index e which 
is 1 for a particle A and —1 for A, so that we 
express the correlation functions as 

oo „ 
L=0 ^ ' m-.-iL'' 

X \K,^,^{e„e,)\^g^ , (4.11) 

where we have introduced the Lx L matrices with 
entries 

K+^{e,(3) = M{e,p), 
K^+{e,p)^M{p,e), (4.12) 
K++{e,p) = K--{e,p) = Q . 

The last equation in 14.1211 ensures that only the 
terms with L = 2n occur in (|4.11(l . Notice that the 
dependence on the statistics in H4.11|l only reduces 
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to taking the permanent or the determinant of the 
same matrix. According to the theory of Fredholm 
integral operators (see e.g. []) the above expression 
is equal to 



G!f^'W=Det(l + 2^K)- 



(4.13) 



so that the correlators of operators having the 
same non-locality phase in the free fermion and 
free boson theories are the inverse of each other. 
This conclusion was first obtained in [] studying 
the deformation theory of differential equations. 
The result (|4.2|l follows from the fact that the 
operators Va{x) have a non- locality phase e^"" 
around the solitons at the sine-Gordon free fermion 
point (see e.g. []), so that Hj{x) and Vj/N{x) have 
the same non-locality with respect to the corre- 
sponding particles. Since j runs between 1 and 
— 1, we are actually working with < a = 
3 IN < 1. 

Once Eq. 1)4. 2() has been obtained by using the 
large distance expansion of the correlators, it is an 
interesting check to see how this inversion relation 
between the two correlators works in the short dis- 
tance limit. Both in the bosonic and fermionic 
case the operators ^a{x) satisfy the OPE 



($„(a;)$„-(0)) - \x\ 



+ (4.14) 



with 



(s) 

Xa being the scaling dimensions. In the bosonic 
case the index a + a' is taken modulo 1. The 
scaling dimensions can be computed through the 
formula [] 

Xi^^ = d^^{Q{x)^c.mconnected, 

(4.16) 

where Q{x) is the trace of the energy-momentum 
tensor. Since the only non-zero form factor of this 
operator in the free theories is 



X^'+X^^-X^^, , (4.15) 



{0\e{0)\A{e)A{l3)) ^ 27rm2 
one easily finds 



-i sinh ■ 



13' 



(4.17) 



in agreement with the conjugation properties fi* ~ 
fiN-j and V* = V-a- These results coincide with 
those of conformal field theory with 'twist' fields 
(see []). It follows 



r 



-2aa' , 



(4.19) 



.(+) 



2aa' , a + a' < 1 

2[aa' + 1- {a + a')] , a + a' > 1 . 

(4.20) 

The agreement for 1 < a -I- a' < 2 is obtained by 
observing that in this range of a -I- a', the leading 
short distance term in the fermionic case is not the 
one in (|4.14|l but the first off-critical contribution 

^^/2 J d^y{Vc.{x)Va.mViiy) + V-iiy)])f,^o , 

(4.21) 

which indeed behaves as +i-{a+a )] .^^jj^g^ 

X 0. 

In general, the two-point functions we are 
dealing with can be expressed as 



where Ta^a'{t) is given by [] 
1 

^ Jt/2 



(4.22) 



(4.23) 



(apx)'-4sinh^X- 



(a - a')' 



tanhx 



with xip) satisfying the differential equation 

dnX+-dnX = 2sinh2Y+.^^^ — tanhr (1-tanh^ 
P 

(4.24) 

with suitable boundary conditions. The short dis- 
tance behaviour for a -I- a' < 1 is given by 



(4.25) 



with an amplitude that can be deduced from the 
work of Ref. [] on vacuum expectation values in 
the sine-Gordon model: 



Cry rv' 2 



(4.26) 



exp < 2 



dt 
~ 
}■ 



sinh at cosh(a -I- a')t sinha'i 
sinh^ t 



Finally, let us discuss the case in which we 
have a Z2 symmetry. This is somehow special 
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since a broken phase with two degenerate vacua 
can be reahsed in terms of a neutralhoson. In fact, 
a disorder operator with non- locality factor —1 is 
present also in the theory of a neutral bosonic free 
particle. It is easy to check that, for S = 1 and 
a = 1/2, the form factors H4.6|l imply the factori- 
sation 



Ml = M(i) X M(2) , 



N 



(4.27) 



where is the disorder operator with scaling 
dimension = 1/8 associated to the neutral 

boson Aj entering the decomposition A = {Ai + 
iA2)/^, i = (Ai ~ iA2)/V2. 

The fermion-boson correspondence observed 
above for charged particles has an analogue in the 
neutral case, and the correlation function 



G 



(+) 

1/2,1/2 



1/2 



G(t) = (/io-)(.T)Mo-)(0)) = 

(4.28) 

can be related to correlators computed in the the- 
ory of a free neutral fermion, i.e. in the Ising field 
theory without magnetic field. To see this, let us 
recall that the correlators of the spin and disor- 
der operators in the unbroken phase of the scaling 
Ising model can be written as [] 



r±(i) = {fi{x)m) ± H^)^m = (4.29) 
= exp 



where the functions x and T are those of Eqs. (|4.23|l . 

(O^ . Hence, it follows from (ji:^ and 

that 

G{t) = [T+{t)T^{t)]-^'\ (4.30) 

Concerning the short distance behaviour of this 
correlator, the power law (|4.25|) acquires in this 
case a logarithmic correction due to the 'resonance' 
with the leading off-critical contribution 14.21|l . 
Such a contribution to the correlator {Va{x)Va{Q)) 
behaves as —Ca,at^^'^[^ + (2 — 4q:) Ini] in the limit 
a 1/2, so that 



limGiy^ 1/2(0 = lmG-2(i) ^ Bt'/' Hl/t) , 

(4.31) 



with 



B = -4ReSa=i/2 C„,a = 0.588353. 



(4.32) 



and this amplitude coincides with that of the prod- 
uct T+{t)T-{t) in the Ising model. 



The considerations of this section have recently 
been extended to the case of massive ghost theo- 
ries [] . In all these cases the correlation functions 
of disorder operators can be computed exactly for 
free theories and are related to the same differen- 
tial equation (|4.24() . In absence of interaction, the 
correlators are non-trivial due to the non-locality 
between the operators and the particles. The re- 
summation of the spectral series for interacting 
theories remains a challanging task even in pres- 
ence of integrability. 

5. Non-integrable Deformations 

The free bosonic and fermionic theories discussed 
above can also be regarded as describing phases of 
spontaneously broken Zn symmetry. In this dual 
vision, the excitations are free kinks \Kj_j±i{9)) 
interpolating between two adjacent vacua |0j) and 
|Oj±i) (the indices are taken modulo N). Let us 
denote by ^k/N, fc = 0, 1, . . . , — 1, the operators 
we are interested in. They correspond to the expo- 
nential operators Vk/N in the fermionic case^, and 
to the operators Ufe, dual to the disorder operators 
Hk, in the bosonic case. These operators create 
multikink excitations with zero topologic charge, 
i.e. starting and ending in the same vacuum state. 
Form factors on kink states were discussed in []. 
Consider for the sake of simplicity the two-kink 
matrix elements 



(01-02) ^ (O,|$,./^(O)|if„±i(0i)i^,±i,,(02)) 



(5.1) 



satisfying the equations 



^;±(0-t-2^7^) = F^^,^,{-e), 
Rese^,^Ff^,{e) = I [(0,|$fc/^|0,) + 

-(Oj±l|<i>fc/Ar|0,±i)] . 

Since the generator of Zn transformations {fl^ 
1) acts on states and operators as (w = e^*'^/^) 

n|i^„-+i(0i)if,+i,,+2(02)...) = 



''In this case the N degenerate vacua can be identified 
with those of the periodic potential in 14.31 identified mod- 
ulo N. 
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the above form factor equations can be rewritten 
as 

F^%{e + 2in) = Sw^''F^,{e), 
Rese=i^F^k{0) = i (1 - a;±'=)(0,|$fe/jv|0,) . 

Once the identifications A < — ifj^j+i, A < — > 
Ki^i-i are made, these relations are equivalent to 
the previous equations for the Form Factors with 
n = 1. This correspondence is easily extended to 
all values of n and leads to the same form factors 
and correlation functions discussed for the unbro- 
ken phase. 

The introduction of a 'magnetic field' pointing 
in the k direction and breaking explicitely the Zn 
symmetry corresponds to adding to the free action 
a term 

h j cfx^k{x), (5.2) 

where 

'i^k{x) = -^u;-''^i/Nix) . (5.3) 

The first order corrections to the energy density 
Sj of the vacuum state \0j) and to the mass mjj±i 
of the kink Kjj±i can be computed according to 
the Form Factor Perturbation Theory proposed in 

D 

5ej ~ /i(0j|*fc|0j) = hvSy,k , 

Smlj^, ~ h{0j\^kmKj,j±i{i7r)Kj±,j{0)) . 

It follows from the residue equation written above 
that 

This formula implies that the correction to the 
mass of the kinks interpolating between the vacua 
|Ofe) and |Ofe±i) is infinite. This divergence sim- 
ply reflects the fact that these kinks become un- 
stable because the magnetic field removes the de- 
generacy of the vacuum jOfe) with the two adja- 
cent vacua. These kinks are then confined in pairs 
Kk,k±iKk±i,k and the confinement gives rise to a 
string of bound states with zero topologic charge, 
as discussed in more details in []. 
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